Single photo-effect on helium-like ions in the non-relativistic region 
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Abstract 

We present a generalization of the pioneering results obtained for single K-shell photoionization of H-like 
ions by M. Stobbe [Ann. Phys. 7 (1930) 661] to the case of the helium isoelectronic sequence. The total 
cross section of the process is calculated, taking into account the correlation corrections to first order of 
the perturbation theory with respect to the electron-electron interaction. Predictions are made for the 
entire non-relativistic energy domain. The phenomenon of dynamical suppression of correlation effects in 
the ionization cross section is discussed. 

PACS numbers: 32.80. Fb, 32.80.-t, 31.25.Eb 
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1. The single ionization of light atomic systems by photon impact is one of the fundamental 
processes, which is being persistently investigated during last decades 

10,Q[j]. Although the non-relativistic problem for single photoeffect on H-like atom was solved 
analytically by M. Stobbe already in 1930 [13J , the further generalization on the case of two-electron 
targets, such as, for example, neutral He atom, is a non-trivial task. Since the strength of the 
electron-electron interaction is comparable to that of the electron-nucleus interaction, the single- 
particle treatment cannot be appropriate. For theoretical calculations of ionization cross sections, it 
is usual to employ sophisticated numerical methods dealing with highly correlated wave functions. 
This allows one to take into account the electron correlations beyond the independent-particle 
approximation . 

In this Letter, we evaluate the dominant contribution of correlation effects to the cross section 
for single K-shell ionization of He-like ions using the non-relativistic perturbation theory. To 
zeroth-order approximation, it is assumed that the atomic nucleus is the external source of the 
Coulomb field, while the interaction of electrons with each other is neglected. The non-interacting 
electrons are described by the Coulomb wave functions for the discrete and continuous spectra 
(Furry picture). The electron-electron interaction is treated within the framework of perturbation 
theory. The latter exhibits fast convergence even for small values of nuclear charge Z [14]. To 
first order of perturbation theory, the problem of single ionization of He-like ions is reduced to 
the evaluation of one-photon exchange diagrams. The correlation corrections to the cross section 
arise due to modification of the binding energy and wave functions of the initial and final states. 
Using this approach, we have already deduced the universal scaling behavior for the ionization cross 
section at high photon energies [3] . The formula obtained previously has a simple analytical form, 
but it is valid only in the asymptotic non-relativistic limit. This Letter reports on the extension 
of results of our work [15| on the entire non-relativistic energy domain. 

2. The non-relativistic problem of single ionization of a K-shell bound electron by photon impact 
involves the following quantities: the momentum k, the energy u = \k\ = k, and the polarization 
vector e of an incident photon, the binding energy / = rf/{2m) and the average momentum 
7} = maZ of a K-shell electron, where m is the electron mass and a is the fine-structure constant 
(h = 1, c = 1). The outgoing electron is characterized by the energy E p = p 2 /(2m) and the 
momentum p at infinity. We shall employ the Coulomb gauge, in which (e • k) = and (e* • e) = 1. 
The parameter aZ is supposed to be sufficiently small, that is, aZ <C 1. 

The amplitude for single photoeffect on a bound K-shell electron is described by the Feynman 
diagram depicted in Fig. [2(a). The energy-conservation law reads E p = u> — I. Accordingly, the 
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non-relativistic photoionization can proceed at the photon energies I ^ u> <C m. Using the one- 
electron Coulomb wave functions for the discrete and continuous spectra, Stobbe has obtained an 
analytical expression for the cross section of single K-shell photoeffect, which can be cast into the 
following form [lij 

+ _ 2 2 9 tt 2 expC-^cot- 1 ^) 

- ««o 3Z 2 (1 + e _ 2)4[1 _ exp(-27r£)] ' { } 

where ao = l/(ma) is the Bohr radius. In Eq. ([I]), we have introduced the dimensionless parameter 
£ = v/P- The quantity £ _1 has the meaning of the momentum p of the outgoing electron, which 
is calibrated in units of the characteristic momentum r/. Due to the energy-conservation law, £ 



is related with the dimensionless energy of the photon e 7 = uj/I according to £ = 1 / y/e^ — 1 or, 
equivalently, £ = where e v = E p /I is the dimensionless energy of the outgoing electron. 

Equation (pQ) is valid in the entire non-relativistic energy domain 1 ^ e 7 <C 2(aZ)~ 2 . In the 
derivation, the dipole approximation has been employed, which implies that the photon momentum 
k is negligible compared with the electron momentum p. 

In the non-relativistic photoeffect, one is often interested in the asymptotic energy domain far 
beyond the threshold, which is characterized by 1 < e T <C 2(aZ)~ 2 . In this case, Eq. (pQ) can be 
expanded with respect to the small parameter £ <C 1. The leading term in the £ expansion is given 
by Q 

4 = ««o^S ?/2 ' (£«!)■ ( 2 ) 

In fact, Eq. ([2]) is obtained using the Born approximation, which describes the wave function of the 
outgoing electron in terms of a plane wave. Since Eq. (pQ) involves also the parameter 7r£, which 
originates from the normalization factor of the Coulomb wave function of the continuous spectrum, 
the convergence of the £ expansion is sufficiently slow. 

Another limit of expression ([T]) corresponds to the ionization threshold (e 7 — ► 1). In this case, 
the cross section reaches the constant value 

2 9 7T 2 

4 = aa o^2 > (£»!)> ( 3 ) 

where e ~ 2.71828 is the Napier-Euler number. 

3. To zeroth-order approximation, the amplitude .A^ -* for the single photoeffect on two-electron 
atomic system in the ground state is described by the Feynman graph depicted in Fig. W[b). In 



3 



the momentum representation, it can be written as follows 

A*> = ^ exp^cot-^) ? (4) 

C = 87T7/iV 7 iV ls iV p (l - iO^p- , (5) 

at9 Ana _ r9 77 3 
AT = iV-r = — 

N 2 = M 

p 1 - exp(-27r£) 

Here £ = 77/p and the dipole approximation (A; = 0) has been again used in the derivation. The 
corresponding expression for the total cross section of the process reads 

a+ = 2a+ , (6) 

where cr^ is given by Eq. (pQ). The factor 2 accounts for the number of electrons in the target. The 
energy-conservation law keeps the same form as for the H-like ion, namely, e p = e 7 — 1. Equation 
§§§ is valid in the entire non-relativistic domain 1 ^ e 7 <C 2(aZ)~ 2 . However, in the case of 
light two-electron atoms, the expression ([6]) is a quite rough approximation, because the electron- 
electron interaction is completely neglected. For example, for the neutral He atom, Eq. ([6]) predicts 
<7q = 3.15 Mb at the ionization threshold (£ — ► 00). The corresponding experimental value turns 
out to be cr+ p = 7.40 Mb at the photon energy to = 24.59 eV [7J]. The significant discrepancy 
between these results is due to correlation effects. Near the threshold, the approximation of non- 
interacting electrons strongly underestimates the experimental cross section. On the contrary, in 
the asymptotic high-energy limit, Eq. © yields an upper bound for the photoionization cross 
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Now we shall take into account the electron-electron interaction using the perturbation theory. 
To first order, the amplitude AS 1 ' for the single K-shell photoeffect on He-like ion is described by 
the gauge invariant set of four Feynman graphs depicted in Fig. [2j The total amplitude for the 
process is given by A — A^ + A^\ where the amplitude (J2|) should be calculated with taking into 
account the correlation correction to the binding energy. More specifically, the energy-conservation 
law now implies e p = e 7 — 1 + 6±, where 5\ = A±/Z is the correction to the binding energy of the 
ground state calculated to first order of the perturbation theory and Ax = 5/4 [17J]. In calculations 
of the amplitude A^ , the correlation correction Si to the binding energy can be omitted, because 
it exceeds the desired level of accuracy. 

Note, that within the asymptotic high-energy domain, the binding energies are assumed to be 
negligibly small with respect to the photon energies. In addition, the wave function of the outgoing 



electron can be approximated by the plane wave, so that the Feynman graph in Fig. E|a) yields the 
dominant contribution to the total cross section for single photoionization, provided the Coulomb 
gauge is employed [lj]]. However, for low energies, all Feynman diagrams depicted in Fig. [2] are 
expected to contribute with the same order of magnitude, while the plane-wave approximation 
becomes inadequate. 

Since the general method for evaluation of the Coulomb matrix elements has been already 
described in more details [3], we focus here on presentation of the explicit expression for the 
amplitude A^ = Yl only. The individual contributions of each diagram in Figs. [2{a)-(d) read 
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Here the function C is defined by Eq. ([5]). In the integral of Eq. ([7|), the prime at the function 
f(x) denotes the derivative with respect to the variable x. We note, that the matrix element 
corresponding to the graph in Fig. [2Ja) involves the reduced Green's function. This requires 
isolation and cancelation of the pole term, which has been done analytically. In Eq. (JSj) , after 



5 



taking the derivative with respect to u, one should set v = 2£, where £ = rj/p. In addition, the 
parameter A of a Yukawa-type screened Coulomb interaction introduced for regularization of the 
integral should tend to zero. In Eq. (|10p . after taking the derivatives with respect to £ and p, one 
should set £ = p = £. Note also that the contribution of the graph in Fig. EJb) describing the final- 
state interaction, contains a logarithmically divergent term ln(A/2p) for A — » 0. However, this term 
drops out in the expression for cross section of the process. This fact reflects a general property 
of the wave function for continuous spectrum, when it is constructed using perturbation theory 
with respect to the Coulomb interaction. More precisely, the logarithmically divergent terms of 
the wave function can be summed up to an overall phase factor exp[i£ ln(A/2p)], which does not 



contribute to the cross section [19|, |2j 
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23]. 



Separating explicitly the contributions due correlation corrections to the binding energy and 
wave functions, the cross section for single K-shell photoionization of He-like ions can be cast into 
the following form 



where crj is given by Eq. Q, £ = 1/y^, e p = e 7 — 1 + 6±, 5\ = A\/Z, and Ai = 5/4. The quantity 
ai(£) = 2Z^Re [A^p /A^^j is a universal function of the parameter £ (see Fig. [3] and Table H]). 

Equation (jlip is valid in entire non-relativistic domain 1 — 5i ^ e 7 <C 2(aZ)~ 2 and represents the 
main result of this Letter. 

The function ai(£) describes the "dynamical correlation", originating from the correlated two- 
electron wave functions both in the initial and final states. Three characteristic domains can be 
distinguished here, namely, the threshold domain (£ 3> 1), the transition domain (£ ~ 1), and 
the high-energy domain (£ ^ 1). As can be seen from Fig. [3l for values £ < 1CP 1 and £ > 10, 
the function ai(£) saturates rapidly approaching the constant limits, while it undergoes significant 
changes within the range 10 _1 < £ < 10. Within the threshold domain, the consistent account of 
the correlation interaction both in the initial and final states becomes especially crucial. In the 

n 

asymptotic high-energy limit, ai(£) tends to the analytical value [15j] 

19 3 

ai = -— + -In 2 ~ -0.6676, (12) 
16 4 

while the correlation correction in the binding energy is negligible for the cross section. In this 
case, the correlation effect results predominantly from the wave functions. In particular, within 
the Coulomb gauge, this is just the contribution from the wave function of the initial state. At 
about £ ~ 0.6036, the function ai(£) tends to zero. Accordingly, the correlation effect in the cross 
section arises due to the binding energy. 
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For some values of the parameter £, the correlation term in the square brackets in Eq. 
vanishes. In this case, the correlation effect due to the binding energy cancels with that due to 
the wave functions. As a result, the cross section fill j) coincides with the Coulomb prediction ([6]). 
In Table HU we consider the particular examples, where the dynamical suppression of correlation 
effects take place. Note, that for neutral helium atom the interpolation of experimental cross 
sections predicted in Ref. [7( yields <7+ = 14.2 kb at the photon energy uj = 308.2 eV. The 
deviation from the corresponding theoretical value = 15.6 kb is exclusively due to higher-order 
correlation corrections neglected in the present consideration. 

In Table IIII) we present a comparison between our predictions according to Eq. (jlip for several 
two-electron targets at different photon energies with the numerical calculations of the work [2] . Bell 
and Kingston used the Hartree-Fock wave functions for the continuum state and many-parametrical 
variational wave functions for the ground state. The results of Ref. [2j are slightly gauge dependent. 
For very light targets, such as neutral He atom, the account for the correlation correction 5\ to 
the binding energy is still not sufficient to reproduce the experimental thresholds. Certainly, one 
needs to go beyond the first-order approximation. However, although consistent calculations of the 
binding e nerg ies have been performed with taking into account of two- and three-photon exchange 
diagrams 



24 



251 ] , the evaluation of ionization cross sections at the same level of accuracy seems to 
be not feasible in the near future. Accordingly, in Table [Hl| we merely employ the experimental 
threshold energies u. As can be seen, the formula (jlip turns out to be in good agreement with the 



numerical results by Bell and Kingston [2j both at the ionization threshold and beyond. 
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TABLE I: For various values of the dimensionless parameter £, the universal quantities ai(£) are tabulated. 
Numbers in parentheses indicate powers of 10. 







ai (0 £ 


ai(0 £ ai(0 C «,(0 e a,(0 £ ar(0 £ 




1.0(- 


-3) 


-0.6676 0.2 


-0.5362 0.7 0.1099 1.2 0.4822 1.7 0.6448 4.0 0.8048 30 


0.8433 


0.5(- 


■2) 


-0.6675 0.3 


-0.4087 0.8 0.2097 1.3 0.5262 1.8 0.6645 5.0 0.8190 50 


0.8438 


1.0(- 


-2) 


-0.6672 0.4 


-0.2693 0.9 0.2955 1.4 0.5634 1.9 0.6816 7.0 0.8313 80 


0.8439 


0.5(- 


-1) 


-0.6580 0.5 


-0.1319 1.0 0.3684 1.5 0.5949 2.0 0.6965 10 0.8378 1.0(2) 0.8439 


0.1 




-0.6309 0.6 


-0.0043 1.1 0.4302 1.6 0.6218 3.0 0.7763 20 0.8424 1.0(3) 0.8440 



TABLE II: For different values of the nuclear charge Z, the parameters £, the dimensionless energies e 7 , 
the Coulomb ionization potentials /, the photon energies u), and the corresponding ionization cross sections 
er + are tabulated. For the considered values of £, the predictions according to Eq. pTjl coincide with those 
according to Eq. ©. 



z 


£ e 7 /(eV) 


to (eV) 


CT+ (kb) Z 


£ e 7 /(eV) 


uj (keV) 


a+ (kb) 


2 


0.4348 5.663 54.42 


308.2 


15.6 


7 


0.4422 5.936 666.7 


3.957 


1.39 


3 


0.4383 5.789 122.4 


708.8 


7.24 


8 


0.4425 5.950 870.8 


5.181 


1.07 


4 


0.4400 5.853 217.7 


1274 


4.15 


9 


0.4428 5.961 1102 


6.569 


0.847 


5 


0.4410 5.891 340.1 


2004 


2.69 


10 0.4430 5.970 1361 


8.122 


0.688 


6 


0.4417 5.917 489.8 


2898 


1.88 


11 


0.4432 5.977 1646 


9.840 


0.570 
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TABLE III: For various two-electron targets, the nuclear charges Z, the photon energies u>, the corresponding 
dimensionless energies e 7 , and the photoionization cross sections a + are tabulated. Our predictions are made 
according to Eq. (fTT|) . The calculations by Bell and Kingston [2| have been performed using different gauges: 
velocity (a) and length (b). Numbers in parentheses indicate powers of 10. 



Target Z uj (eV) e 1 a+ (Mb) w (keV) e 7 a+ (b) 

Eq. (dU) Ref. [2] a Ref. \2] h Eq. $TT$ Ref. [2] a Ref. [2] b 



He 


2 


24.59 


0.452 7.19 


7.15 


7.32 


0.456 


8.375 4.68(3) 


4.40(3) 


4.80(3) 






61.23 


1.125 1.21 


1.07 


1.14 


0.891 


16.37 553 


527 


558 






183.7 


3.375 7.19(- 


^2) 6.64(- 


-2) 7.34(- 


-2) 1.327 


24.37 151 


146 


151 






238.1 


4.375 3.39(- 


-2) 3.15(- 


-2) 3.47(- 


-2) 1.762 


32.37 59.5 


57.4 


59.0 






346.9 


6.375 1.09(- 


-2) 1.02(- 


-2) 1.12(- 


-2) 3.503 


64.37 6.05 


5.74 


6.12 


Li+ 


3 


75.64 


0.618 2.49 


2.57 


2.61 


0.289 


2.361 9.88(4) 


9.56(4) 


9.96(4) 






112.2 


0.917 1.07 


1.02 


1.05 


0.398 


3.250 4.01(4) 


3.88(4) 


4.06(4) 






153.1 


1.250 5.19(- 


-1) 4.98(- 


-1) 5.11(- 


-1) 1.378 


11.25 918 


897 


928 






180.3 


1.472 3.45(- 


-1) 3.33(- 


-1) 3.43(- 


-1) 1.813 


14.81 383 


375 


386 






234.7 


1.917 1.74(- 


-1) 1.68(- 


-1) 1.74(- 


4) 3.554 


29.03 43.2 


41.7 


44.0 


B 3+ 


5 


259.3 


0.762 0.730 


0.743 


0.748 


0.473 


1.390 1.67(5) 


1.65(5) 


1.67(5) 






295.9 


0.870 0.537 


0.528 


0.533 


0.582 


1.710 9.65(4) 


9.53(4) 


9.64(4) 






336.7 


0.990 0.393 


0.388 


0.391 


1.126 


3.310 1.50(4) 


1.48(4) 


1.51(4) 






363.9 


1.070 0.325 


0.321 


0.323 


1.561 


4.590 5.72(3) 


5.64(3) 


5.74(3) 






418.4 


1.230 0.229 


0.226 


0.228 


3.738 


10.99 388 


384 


399 


Q6+ 


8 


739.0 


0.849 0.251 


0.253 


0.254 


3.782 


4.344 2.77(3) 


2.76(3) 


2.78(3) 






1497 


1.719 3.99(- 


-2) 3.98(- 


^2) 4.00(- 


-2) 5.306 


6.094 994 


990 


996 






1878 


2.156 2.13(- 


-2) 2.12(- 


-2) 2.13(- 


-2) 6.830 


7.844 457 


455 


458 






2259 


2.594 1.26(- 


-2) 1.25(- 


-2) 1.26(- 


-2) 9.878 


11.34 144 


144 


144 






3020 


3.469 5.40(- 


-3) 5.38(- 


-3) 5.41(- 


-3) 12.92 


14.84 61.4 


61.1 


61.4 


Ne 8 + 


10 


1195 


0.878 0.153 


0.154 


0.155 


7.068 


5.195 1.05(3) 


1.05(3) 


1.05(3) 






1313 


0.965 0.121 


0.121 


0.121 


9.027 


6.635 498 


497 


498 






2660 


1.955 1.83(- 


-2) 1.83(- 


-2) 1.83(- 


-2) 12.95 


9.515 163 


162 


163 






4129 


3.035 5.20(- 


-3) 5.19(- 


-3) 5.21(- 


-3) 16.86 


12.39 70.9 


70.7 


70.9 






5109 


3.755 2.78(- 


-3) 2.77(- 


-3) 2.79(- 


-3) 24.46 


17.97 21.7 


21.4 


21.5 
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J. 



Ipp — < ^ «— 01, 

(a) (b) 

FIG. 1: Feynman diagrams for the single ionization of a K-shell electron by a single photon. Diagram (a), 
hydrogen-like ion; diagram (b), helium-like ion without taking into account the electron-electron interaction. 



lp p < $ T-^<— ^xs A — * — T *l>u 

i i 
i i 

^Xs *l>ls ^Xs — *— * <— V>15 



(a) (b) 
A: 



^xs < ** T—<— ^xs ^xs — *— t ^1, 

I I 
I I 

ll> p — * *— «— 1p ls 1/> P «— * «— -01, 

(c) (d) 

FIG. 2: Feynman diagrams for the single K-shcll photoionization of He-like ion. Diagrams (a) and (c) take 
into account the electron-electron interaction in the initial state, while diagrams (b) and (d) account for it in 
the final state. The individual contributions of each diagram are gauge dependent, while total contribution 
of all diagrams is gauge independent. 
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FIG. 3: Different contributions to the universal function ai(£) calculated in Coulomb gauge. Dash-double 
dotted line, contribution due to the diagram in Fig. Ufa); dash-dotted line, contribution due to the diagram 
in Fig. EJb); dotted line, contribution due to the diagram in Fig. EJ^c); dashed line, contribution due to the 
diagram in Fig. (5fd); solid line, total contribution of all diagrams. 
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